The generalised Dirichlet to Neumann map for moving 
initial-boundary value problems 

A.S. Fokas* and B. Pelloni** 



*Department of Applied Mathematics and Theoretical Physics 
Cambridge University 
Cambridge CBS OWA, UK. 
t.fokas@damtp. cam. ac.uk 



**Department of Mathematics 
University of Reading 
Reading RG6 6AX, UK 
b.pellom@rdg. ac. uk 



February 6, 2008 



Abstract 



We present an algorithm for characterising the generalised Dirichlet to Neumann 
map for moving initial-boundary value problems. This algorithm is derived by combin- 
ing the so-called global relation, which couples the initial and boundary values of the 
problem, with a new method for inverting certain one-dimensional integrals. This new 
method is based on the spectral analysis of an associated ODE and on the use of the 
d-bar formalism. As an illustration, the Neumann boundary value for the linearised 
Schrodinger equation is determined in terms of the Dirichlet boundary value and of the 
initial condition. 

1 Introduction 

We present a methodology for characterising the generahsed Dirichlet to Neumann map for 
linear evolution PDEs posed on domains whose boundary varies with time. Consider, as an 
example, the following domain in the {x,t) plane, see figure [T] 



where T is a positive constant and l{t) is a given monotonic, twice differentiable function. 
For economy of presentation we assume that l{t) satifies the following: 



D : < i < T, l{t) <x <oo, 



(1.1) 



l{t) e C^[0,T], l"{t)>0, l{0)^0. 



Let the scalar complex- valued function q{x,t) satisfy the boundary value problem 



iqt + qxx = 0, {x,t) e D, 

q{x,0)^qo{x), 0<x<^; = /o(i), 0<t<T. 



(1.2) 
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Figure 1: The domain D in the [x, t) plane 



We assume that the initial condition qaix) is a sufficiently smooth function, decaying as 
X —^ OO, that the boundary condition fo{t) is sufficiently smooth, and that these two 
functions are compatible at the origin, qo{0) = /o(0). 

There exist several different integral representations for the solution of initial-boundary value 
problems such as those defined by p.2p . These include the representation constructed by 
using the classical Fourier transform in the x variable, as well as the novel representations 
presented in [6l[7]. However, these representations are not effective, because they involve 
unknown boundary values. For example, for the initial-boundary value problem ()1.2p these 
representations involve the unknown function qx{l{t),t). 

Thus, the key issue is the characterisation of the generalised Dirichlet to Neumann map, 
namely the characterisation of the unknown boundary values in terms of the given initial and 
boundary conditions. For the initial-boundary value problem (jl.2l) this means computing 
qx{l{t),t) in terms of {(7o(a;), /o(i)}- 

Here we present an algorithm for constructing the generalised Dirichlet to Neumann map. 
For the initial boundary value problem (|1.2p this algorithm yields the unknown function 
qx{l{t),t) = fi{t) through the solution of the Volterra integral equation 



/iW = |^^ J(s,i)/i(s)ds+ 



(1.3) 



+ 



30F 



t/4 



e •*('-=■) 1 f i(l{t)-a,)^ 

fo{s)ds + — / e 4t q^{x)dx 

Vt- s Vt Jo 
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for < t < T, where the kernel J(s, t) is defined by 

2(s-t) [Jl'(s)/2 Jo 

(1.4) 

and 

G{s,t) = e'^(«-*)-^«^)-'(*)). (1.5) 



We note the both integrals on the right hand side of equation (|1.4p are well defined. In 
particular, the second integral involves an exponential which decays as k ^ oo. Indeed, the 
real part of the exponent is 

-k{s - t) (^l'{s) - '^^r^^Y^^ = -Hs - t){l'{s) " l'{cr)), s<a<t. 

Recalling that l"{t) > 0, it follows that l'{s) — I' (a) < 0, and since k > and s — t < 0, the 
real part of the exponent is negative. 
The algorithm involves three steps: 

1. Assuming that the solution exists, derive the so-called global relation, namely the relation 
that couples the initial condition with all boundary values. This step is elementary and it 
involves only writing the given PDE in a proper divergence form, and applying Green's 
theorem. For example, equation (jl.2p can be rewritten in the form 



6-*'^'^+*'= *g(a;, i) - 6"*'^'^+*'= *(ig^(x, t) - fcg(x, t) =0, k e C. (1.6) 

J i L J X 

Then an application of Green's theorem in the domain D yields 

e'''''-'^^(-'Ux{l{s),s)ds^ / e'^''-''''^-'\k-l'{s))fo{s)ds 
Jo 

POO 

+ qo{k) - e'''''^ e~'''''q{x,T)dx, Im(A:) < 0. (1.7) 

Jl(T) 

2. Consider the integral involving the unknown boundary values and derive a general formula 
for the inversion of this type of integrals. This step involves the spectral analysis of an 
appropriate ODE. For example, for the initial boundary value problem (|1.2p . this ODE is 

Ht{t,k) + {ik^ -ikl'{t))ii{t,k)^kf{t), 0<t<T, keC. (1.8) 

Using this ODE it is shown in Section 2 that if F{k) is defined in terms of f{s) by the 
integral 

rT 

F{k)= e'^''-'^^'^'^f{s)ds, fceC, (1.9) 



then f{t) can be obtained in terms of F{k) through the solution of the following Volterra 
integral equation 



2 f n rt 



fit) - — / e-^'^ '+''''^'^F{k)kdk + — / J(s, t)f{s)ds, 0<t<T, (1.10) 
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where J{s,t) is defined by equation (|1.4p and 9f22 (t) is given by 

dfl- (t) ^ {k e R : k<l'{t)/2}U{k^kii + iki : k^ ^ l'{t)/2, ki < 0} (1.11) 
with the orientation shown in Figure [21 



Figure 2: The curve (Jfij (t) in the fc plane 



3. Use the global relation |_?. 7[ ) anrf t/ie inversion formula hl.lOfl to derive a Volterra integral 
equation for the unknown boundary values. The global relation (|1.7p is of the form (jl.9p 
where /(t) and F(fc) are replaced respectively by ifi{t) and by 

qo{k)-e'^^^ e-'''''q{x,T)dx+ e*^'^-*'='(^)(fc - Z'(s))/o(s)ds. (1.12) 

Jl(T) Jo 

Furthermore, the global relation is valid for Im{k) < 0, therefore it is valid for k e dft2{t). 
It is shown in Section [3] that equation (|1.3p follows from equation (jl.lOp by replacing in 
equation (|1.9p /(s) and -F(fc) by i/i and by the expression in (|1.12p . 

It turns out that the second term in (|1.12p yields a zero contribution, which is consistent 
with the evolutionary nature of equation (jl.2p {q{x, t) cannot depend on the future time T). 
Although the generalised Dirichlct to Neumann map is obtained under the assumption of the 
existence of a unique solution, this map can be justified a posteriori without this assumption, 
see Theorem 1.1 of ^7,. 



2 The spectral analysis of ordinary differential equa- 
tions and the inversion of complex integrals 

Proposition 2.1 Let F{k) be defined in terms of f{t) by equation il.9\) . where f{t) is a 
sufficiently smooth function. Then f[t) satisfies the Volterra integral equation il.lO\) . 

Proof: (a) We first treat the ODE (II. 8p as an equation which defines fi{t, k) in terms of 
f{t) and we seek a solution which is bounded for all values of the complex parameter k. 
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By integrating with respect to t from either or T we find the following two particular 
solutions of (|1.8p : 

fxi{t,k) - re^('='(^-*)-'=('(^)-'W)fc/(s)ds, (2.1) 
Jo 

lJL2{t,k) = -J e'('='(''-*)-'=('('')-'W)fc/(s)ds. (2.2) 

The functions /ii and /i2 are entire functions of fc, which are bounded, respectively, in the 
domains Qi and Q2 defined by 

n^{t) = {keC: ki>0, kR<'^}U{keC: ki <0, kji>'-^, 0<t<T}. 

(2.3) 

l'(T) I'lt) 
n2{t) = {keC: ki>0, kR> ^^}U{keC: ki <0, kR<^, 0<t<T}. 

(2.4) 

The domains ^i{t) and f22(i), which are depicted in figure[3l are determined by the real part 
of the exponent of the exponential term appearing in equations (|2.ip . (|2.2p . which equals 

g-(s-t)fci(2fcH- "°]:'/*' ) ^ ^_(s_t)fc,(2fe„-;'(r))^ 

where r is in the interval bounded by s and t. 
For fii, s — t < 0, thus fii is bounded if and only if 

ki{2kR-l\T)) <0, 

i.e. 

1'(t) 1'(t) 
{ki >0,kn< -^} or {kj <0,kn> 1^}, 

for every r in the interval < s < t < t. Taking into account that l'(t) is an increasing 
function, the above inequalities yield the definition of fli. Similarly for fi2 and ^2it)- 
A solution of the ODE (|1.8p which is bounded in the domain 

I'iQ) l'(T) 
ns^{keC: ki>0, ^<kR< AJ}, (2.5) 

is given by 

nt 

H3it,k,k)= e*('^-'(^-*)-'=('(^)-'W)fc/(s)ds, (2.6) 

where the function S{kji) is defined on the real interval [^-^, ^-^^1 by 

V{0) V{T) l'{s) 
'^(fcfl) : [-^, — [0,rj, S{kR)^s^=^kR^^^. (2.7) 

In order to prove that the function is bounded in l^a, we distinguish two cases: 
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Figure 3: The domains Qi, Q2 and in the fc-plane and the contours ri2,ri3 and T23 



• < S{kR) < t 

In this case, s — t < and since kj > 0, we need to prove that kj^ < /'(t)/2, which 
follows from the fact that Z'(r)/2 > /'(s)/2 > fc^, 

• t < S{kii) < T 

In this case, s — t > and since fc/ > 0, we need to prove that ku > 1'{t)/2, which 
follows from the fact that 1'{t)/2 < l'{s)/2 < kR. 

We emphasise that the function fj,3{t,k,k), in contrast with the functions ^i{t,k) and 
^2{t, k), involves fc^,, hence this function dependes on both k and k. 

Integration by parts of the equations defining the /ij's implies the following asymptotic 
behaviour 

1" 



fi,=Oi^-j, ken,, fc^«3, j = l,2,3. (2.8) 

Equations (|2.ip . (|2.2p and (|2.6p define a function fi{t, fc, k) in terms of f{t). 

(b) Using the fact that the function /i(t, fc, k) is bounded in the entire complex k plane, 
including infinity (see equation (|2.8p ). it is possible to find an alternative representation for 
this function using the Pompeiu ( also known as Cauchy-Green, or d-bar) formula [1], 

,it,k,-k)^^f ^^^dX + ^ff Mt,\X)dXAdX^ 0<t<T,keC, 



2m 7^.(4) A - fc 2m J 7o(t) dX A - fc ' 

(2.9) 

where dX A dX = -~2idXiidXi , ^{t) denotes the contour along which the function ji has a 
"jump" discontinuity and Q,{t) is the domain where dfi/dk 7^ 0. Before computing the 
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Ml - = 



M2 - 



relevant jumps, we note that 113 coincide with /ii on the half-hne {kj > 0, kji — Z'(0)/2} 
and with ^2 on the half-hne {kj > 0, kj^ ~ l'{T)/2}. This is consequence of the definition 
of S{kii), which imphes that S{l'{s)/2) = s for all < s < T. Thus, the relevant jumps are 
given by the expressions 

= / A{k,s,t,)ds, fc e ri2, 

/ A{k, s, t, )ds, k e ri3, 
Jo 

j A{k,s,t,)ds, k e r23, (2.10) 

where we use the following notation (see figure [3]): 

A{k,.s,t) = e^k\s~t)-^k{^s)-l(ty)f^J^f^^■^^ (2.11) 

Tuit) = {fc7 = 0, /cflG (-(X3,^)}U{fcfl,= ^, fc/G (0,-cx))}U 

U{fc/ = 0, fci^ e (oo,Z'(r)/2)}, 
l'(T] l'(t) 

ri3(t) = {fc/-o, (A^,^)}, (2.12) 

T2,(t) = {fc.^O, kae (^,^)}. 

The direction of integration is depicted in figure[31 and it is also indicated in equations (|2.12p : 
for example, k^ g (— 00, Z'(0)/2) indicates that the integration is from —00 to Z'(0)/2. 
Equation (12. 6|) implies 



a^^aM^ = -^He-^(^-*)-('(^)-'(*»/(.)}...(..) = ~'-^A[k,S{k,),t). 

(2.13) 

Hence, using equations (|2.10p and (|2.13p in equation (|2.9p . we find the following alternative 
representation of /i: 

^(i, /c, fc) = / (/ ^(A,s,t)ds ) -^ + -1^ / [ [ ^ ''\{\s,t)ds '^^ 



27ri 7ri2 l^'o I X — k 2ni 7ri3 V-'o I X — k 

1 f ( f^'^'"'^ .,\dX Iff dS(XB) r.,. N ,dXhdX,^ 



(c) The representation of the function ^{t,k,k) defined by equations (|2.ip . p.2p . (|2.6p 
involves f{t), while the represnetation defined by equation (|2.14p involves various integrals 
of A{k, s, t). Thus there exists a relation between f{t) and these integrals of A. The simplest 
way to obtain this relation is to consider the large fc asymptotic behaviour of /i(t, fc,fc). 
Equation (|2.14p implies 

^{t, fc, fc) = + o(y], Mo - lim (fc^). 



fc / fc^oo 
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Substituting this expression in the ODE (|1.8p we find = f{t), thus 



■S(kii) 

A{k,s,t)ds ] dk ~ / ( / A{k, s,t)ds ] dk 
ri3 



^ ''\{k,s,t)ds^ dk + J J ^^^BlA{k,S{kR),t)dkAdk^. (2.15) 
Equation (|2.15p was first derived in [7] (see equation (3.2) of [7]). 

(d) We will now show that equation (|2.15p can be transformed into a Volterra integral 
equation. For this purpose, we split the domain fl^ in the form fl^ — fl'i^\t) U fl'"^\t) where 

ni^\t) = {fc e : l'{Q)/2 <kR< l'{t)/2}, n'^^\t) = {fc e : l'{t)/2 < kn < l'{T)/2}. 

(2.16) 

(2) 

In the domain ilg ' (t) we use the complex form of Green's theorem, which states that 

a/i3(fc,fc)^^^^- ^ r ^^^k,k)dk, (2.17) 



o("'(t) dk Jdfif\t) 

f2) ('21 

where the boundary Qfig ' (t) of Jig ' (t) has counterclockwise orientation. Recalling that 
S{kfi) = t when 2^/? = Z'(t), it follows that the contribution of the half line {kji = 
l'{t)/2, kj > 0} to the right hand side of equation (|2.17p vanishes, hence 

^^^^^^ A{k, Sikn), t)dk Adk^- f f ^^^4r^dk Adk = 



f2(^'(t) J -'f^f (t) 

rl'iT)/2 . 
= / / A(fc,s,i)dsdfc+ / / A{k,s,t)dsdk, (2.18) 

Jl'{t)/2 Jsikn) Ji(T)Jt 

where the half line 1{T) is defined by 

1{T) = {k:kR^ l'{T)/2, 0<ki < oo}. (2.19) 

The integrands of the line integrals appearing on the right hand side of equation (|2.15p are 
shown in figure 01 while the integrands of the line integrals appearing on the right hand side 
of equation (|2.18p are shown in figure [5l 

Adding the integrands in the respective line integrals, we find the contributions depicted in 
figure (O. 

The above analysis, together with the identity 

dSikn^^dSp)^ dfc A A = -2*dfc«dfc,. 
dk 2dkji 
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Figure 4: The integrands in equation (|2.15p 




Figure 5: The integrands in equation (|2.18p 



impHes that equation (|2.15p can be rewritten in the following form 

2nf{t) = / (/ A{k,s,t)ds] dk+ [ if A{k,s,t)ds] dk 
Jdfi-(t) \Jo J Jdn+(t) \Jt J 

l'(t)/2 I fS(kR) \ { (°° ( \ 

\ A(k,s,t)ds\dk^ I / / A(k , s ,t)ds \ dk 

;'(0)/2 \J0 I \Jl'{t)/2 \Jo J 



' fdSj^ I A(fc,^(fcfl),i)&« I dki, 



V(0)/2 



dkn 



(2.20) 



As a consequence of splitting the integrals originally appearing together on the right hand 
side of equation (|2.14p . the resulting individual integrals in equation (|2.20p have a singularity 
as fc — > oo. Hence the last two integrals on the right hand side of equation (|2.20p must be 
considered together (this is indicated by the curly brackets around these terms). 
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Figure 6: The integrands in the expression obtained by adding (|2.15[) and (|2.18p together. 
The double hnes indicate that the contribution is spht into two integrals. 



The first term on the right hand side of this expression can be computed in terms of the 
given data F{k), defined in (|1.9p . hence this term is known. 

The second term equals (7r/2)/(t). Indeed, the integrand in this term is bounded and 
analytic in ri^(t), and it behaves like if(t)/k as fc ^ oo. Thus this term equals 

J-k/2 ^ 

The third integral in (|2.20p can be written in a "Volterra form" by exchanging the order of 
integration. Indeed, recalling the definition of S{kfj) and changing variable to cr = S{kjj) so 
that k{= kji) = l'{a)/2, we find that this integral is equal to 

i'(t)/2 / .Sikn) \ rt / nl'{t)/2 \ 

/ A{k,s,t)ds\ dk ^ / / A{k, s,t)dk ] ds 

l'(0)/2 \Jo J Jo \Jl'(s)/2 J 

In the double integral we change the variable kji to s — S{kji) and rename fc/ as k. This 
yields 

^'-''^'^ <^S^tRLA{k,S{kR),t)dkRdki= [ ([ A{^-^+tk,s,t)ds)dk. 
lO Jl'(0)/2 "^fl Jo \Jo 2 J 

Hence equation (|2.20p yields 



-m-f ( 

2 Jdn-{t) \ 



A{k,s,t)ds\dk+ i i A{k, s,t)dk \ ds+ (2.21) 

J Jo \Jl'{s)/2 J 



A{k,s,t)ds] dk-i [ ([ Ai'-^ + ik, s,t)ds) dk\ . 
l'{t)/2 \Jo J Jo \Jo 2 / J 

Setting 

Aik, s, t) = E{k, s, t)fc/(s), E{k, s, t) = Q^k\s-i)-^k{l(s)-l(t)) ^ (2.22) 



and using 



+ ik, s, t) = E{ik, s, t)e-'=("-*)''(^)G(s, t), 
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where G{s,t) is given by (|1.5p . equation (|2.2ip becomes 



^ Jd 



dn~{t) \Jo 



A{k,s,t)ds dk + 



t ( /•/'(t)/2 



\Jl'{s)/2 



E{k,s,t)kdk\ f{s)ds+ (2.23) 



f E{k,s,t)fis)ds\dk-i f (f £;(i/c,s,i)(fc-i^)G(s,i)e-*=("-*)''(')/(s)rfs 
^ /;'(t)/2 \Jo ) Jo \Jo 2 

In what foUows we will rewrite the terms in the above bracket in a Volterra form. In this 
regard we note that this bracket is well defined. Indeed, integration by parts of the inner 
integrals yields 



l'{t)/2 



m 

ik 



dk — i 



m 
k 



dk. 



which shows that the contribution of the singular term 1/fc at infinity cancels out. 
We start by rewriting the bracket as 



Ms,t) + G{s,t)J2{s,t) ~ -G{s,t)l'{s)J2{s,t) 



fis)ds 



(2.24) 



where 



Ji(s,i) 


POO 




e^0+ Jl'{t)/2 




J2is,t) 


poo 

= lim / Q-^i'''i^-t+i'^)+KKs)-m)-k{s-t)i'{s)i^^i^ 




-^-0+ io 




J2{s,t) 


= / E{ik,s,t)e'^^^'-'^^'^^'Uk. 
Jo 


(2.25) 



The kernel in the expression (|2.24p equals 

g(g,t) l{s)-lit) 
2i{s-t) 2{s~t) 

where G{s,t) is given by (|1.5p and 



Ji{s,t) - iG{s,t)J2{s,t) 



(2.26) 



if(s,i) = G(s,i)-e^^("-*)-*^('(")-'(*», Ji(s,i) = 



l'(t)/2 



E{k,s,t)dk. (2.27) 



Indeed, integrating by parts the expressions for Ji and J2 we find 



Jiis,t) 

J2{s,t) 



2i{s~'t + i0) 



2(s - t) Jl'{t)/2 



E{k,s,t)dk, (2.28) 



2i{s -t + iO) 



lis) -lit) I'is) 



2i(s - t) 2i 



Ms^t). 



(2.29) 
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Substituting the above expressions in (|2.24p we find that the kernel equals the expression in 

Using integration by parts of the second term on the right hand side of (j2.23p . we find 

H(s t) l(s) - l(t) 
Eik,s,t)kdk = --^^ + y / E{k,s,t)dk. (2.30) 

l'(s)/2 2l{s-t) 2{s-t) Ji,(^s)/2 

Adding the expressions in (|2.26p and (|2.30p we find that the equation (|2.23p can be written 
in the form 

^/W- / (/ A(k,s,t)ds]dk+ f J{s,t)f{s)ds, 

^ Jdil^it) \Jo J Jo 

where J(s, t) is the kernel given by (|1.4[) . 



Remark 2.1 The direct evaluation of the first integral appearing in expression (|1.4p . by 
completing the square in the exponent, yields 



oc 



l'(s)/2 Jl'(s)/2 



Setting A = [fc - ^-m],/t—s, and \o{s,t) = [?'(s)/2 - 

written as 

(i(t)-Ut))^ 



this can be 



OO 



E(k R f)dk^ ^ I e-'^ dX. (2.31) 

l'(s)/2 Vt-S J\o{s,t) 

The integral on the right hand side of equation (|2.3ip can be expressed in terms of the 
Gamma function, and of the imaginary error function Erfi(z). Indeed, 

e-'^'dX= f e-'^'dX+ I e-'^'dX 



i + 1 [tt 



2 V 2 ' "V > // 2 



3 The Dirichlet to Neumann map for the hnear Schrodinger 
equation 

Proposition 3.1 Let the complex-valued scalar function q{x,t) satusfy the following initial- 
boundary value problem: 

iqt + qxx = 0, l'{t) <x <oo, Q<t <T, 

q{x, 0) = qoix), l'{t) < x < OO, , , 

qm),t)^ fo{t), 0<t<T, ^''■'> 
90(0) = /o(0) 

where l{t) G C^[0,r], l"{t) > 0, l{0) = 0. The Dirichlet to Neumann map for this 

problem is characterised by the linear Volterra integral equation il.3\) . 
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Proof In order to obtain the linear integral equation satisfied by fi{t) = qx{l{t),t) we 
must replace in equation (|1.10p the function f{t) by ifi{t), and the function F{k) by the 
expression in (|1.12p . Multiplying the latter expression by k we find three terms. 
The first term in (|1.12p . multiplied by k, equals 

) /"OO 

e-'''''qo{x)dx^-iqo{Q)-i e-'^"" q'o{x)dx. (3.2) 
The second term multiplied by j^Q-'^k^^+'^'^iit) jg 

e-^K--iiT))q^^^T)dx. (3.3) 

The exponential multiplying the above integral is bounded in the domain fij^ (t) , while the 
integral in ()3.3|) is bounded and analytic for Im(fc) < and is of order O (^) as /c — > oo. 
Thus, the application of Jordan's lemma (after a suitable change of variables) implies that 
the integral of (|3.3p along dfl2{t) vanishes. 
The third term in (|1.12|) . multiplied by k, equals 

e^'='^-*"(")(fc2-fcr(s))/o(s)ds = -ie^'='^-^'='(^)/o(T) + i/o(0) + i / e^'='^-*^-'(")/^(s)ds. 
Jo 

(3.4) 

Adding the terms on the right hand side of equations (|3.2[) and p.4[) (using the compatibility 
go(0) — /o(0)), multiplying the resulting expression by e~*'^ t+iki{t) integrating it with 
respect to k along the contour 917^ (t), we find 



ldn-(t) 



37r 



t 



Mds- / e-'''-q'^{x)dx 



dk 



+ |- / J{s,t)f,{s)ds, 0<t<T. (3.5) 

Indeed, the integral involving Jq{T) as well as the integral involving J^"^ e*'^^^~''^''^*^/Q(s)(is 
vanish, because the term exp[— ifc'^(t — s) + ik(l{t) — l{s))] is bounded and analytic in ^2^^) 
whenever s — t > 0. 

After exchanging order of integration, the integral over k on the right hand side of equation 
(|3.5p can be computed explicitly, and when this is done equation (13. 5|) becomes equation 
(fL3| . Indeed, 

— I e'fc^(^-*)-*'=W^)-'(*))dfc = f — / e-^^'dA^ e^2^^ (3.6) 

and 

(\ ■ (i(t)-x)^ 

A / e-'^'dA ^' ; , (3.7) 
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a(t, s) 



b{t,x) 



Figure 7: (7a) The domain fij (^b) The domain {t, x) 



where the domains Q2 (^j ^^'^ ^2 •'^) ^^'^ depicted in figure ([7]) and the fimctions a{t, s) 
and 6(t, x) appearing in the figures are given by 



a{t, s) 



l'{t) 



m - m 



t - s 



b{t,x) 



l\t) 



l{t) - X 



Since t > s it foUows that ^-^zj^ = I'ir), s < t < t, thus l'{t) - 1'{t) > and hence 



a(t, s) > 0. Similarly, 



nt)-^-^^ = nt)-nr) + ^, r<t, 



thus b{t, x) is also positive. 

The exponential e~ is bounded in the fourth quadrant of the complex A plane, so that 
both the lines a + ikj and b + ikj, kj < 0, can be deformed to the real axis. Hence both the 
A integrals equal 



2 ,2 2 



Stt 



2 / e-*-^ dX 



-i7r/4 



QED 



4 Conclusions 

The results presented in this paper are perhaps interesting in a broader context, as they 
illustrate the implementation of a new technique, which allows one to invert complicated 
integrals, such as the integral defined by equation (jl.9p . This integral is a simple variant of 
the elementary integral 

T{k) = / e^^''f{s)ds, keC. (4.1) 
Jq 

This integral can be inverted by a straightforward application of the inverse Fourier trans- 
form (after a suitable change of variables) 

f{t) = - / e-'^^'^kT{k)dk, 0<t<T, (4.2) 
J ax 
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where dX denotes the boundary of the first quadrant I of the complex fc-plane, with coun- 
terclockwise orientation. It is interesting that small variations of elementary integrals such 
as the integral (|1.9p , apparently have not been investigated until now. Perhaps this is due to 
the fact that the analysis of the integral l|1.9p . in contrast to the analysis of the integral (|4.ip . 
involves functions that are not analytic. Indeed, besides using the Fourier transform, there 
exists an alternative approoach for inverting (|4.ip . which is based on the spectral analysis 
of the ODE 

^j.t{t,k) +ik^fi{t,k) ^ kf{t), 0<t<T, keC. (4.3) 

The crucial difference between this equation and the ODE p.Sp (which is associated with 
the integral p.9p ) is that whereas there exists a solution n{t,k) of equation (|4.3p which is 
sectionally analytic, the solution of (|1.8p involves Hsit, fc, k) for which 8^3/ dk ^ 0. 
The method used in this paper for inverting the integral (|1.9I) has its origin in the paper [4] , 
where it was emphasised that techniques developed for the solution of integrable nonlinear 
PDEs provide a new method for constructing integral transforms pairs. In particular, it was 
shown in that the spectral analysis of the ODE 

fix{x, k) — ik^{x, k) — q{x), a: G M, G C, 

yields the Fourier transform pair. The first nontrivial application of this method appeared 
in the work of R. Novikov [lOj , who was able to invert the attenuated Radon transform by 
using a simple extension of a novel derivation of the Radon transform obtained in [5] using 
the methodology of [1]. It appears that our work, which yields the inversion of a large class 
fo integrals, presents a second nontrivial application of the method of [1]. These integrals 
are precisely the ones that characterise the Dirichlct to Neumann map for moving initial- 
boundary value problems. Given the simple form of these integrals, it is natural to expect 
that they may appear in other applications. The inversion of of the integral characterising 
the Dirichlct to Neumann map for the heat equation is presented in [3]. 
The spectral analysis of equation ()1.9p was first carried out in [6l [7j where the particular 
solutions /xi, fj,2 and /X3 (see equations (|2.1I) . (j2.2p . (|2.6p ) were introduced. However, in 
[6l [7] the inversion formula for f(t) was left in terms of a two-dimensional integral, and 
thus it did not provide an effective way of constructing f{t). The crucial new development 
presented here is the understanding that the double integral can be expressed in terms of 
single integrals, which in turn yield a Volterra integral equation for f{t). 
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